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ABSTRACT 



Recent work has shown that it is possible to systematize quasar spectral diversity 
in a parameter space called "Eigenvector 1" (El). We present median AGN spec- 
tra for fixed regions of the El (optical) parameter space (FWHM(H/?) vs. equiva- 
lent width ratio R FcI i=W(Fell A4570)/W(H/3). Comparison of the median spectra for 
different regions show considerable differences. We suggest that an El-driven ap- 
proach to median/average spectra emphasizes significant differences between AGN, 
and offers more insights into AGN physics and dynamics than a single population me- 
dian/average derived from a large and heterogeneous sample of sources. We find that 
the H/3 broad component line profile changes in shape along the El sequence both 
in average centroid shift and asymmetry. While objects with FWHM(H/3bc)^ 4000 
km s' 1 are well fitted by a Lorentz function, AGN with FWHM(H/3 B c)> 4000 km s _1 
are well fitted if two broad line components are used: a broad (the "classical" broad 
line component) and a very broad/redshifted component. 

Subject headings: quasars: emission lines - quasars: general - galaxies: active 



1. Introduction 

Studies of broad emission line spectra for AGN provide both the strongest constraint on mod- 
els for the nebular physics and kinematic models of the broad line emitting clouds. Ideally one 
would like to have both high signal-to-noise ratio (S/N) and moderate resolution ( ~ 5 A) measures 
for the strongest high and low ionization lines as well as measures across as wide a wavelength 
range as possible in order to better characterize the continuum shape. The former measures have 
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led to the Eigenvector 1 parameter space concept (Sulentic et al. 2000a,b) that is, in part, built upon 
foundations laid almost ten years ago (Boroson & Green 1992; hereafter BG92). The above goals 
are very observing time intensive but recent AGN surveys have provided an alternate method for 
examining quasar spectra from 1000 - 10000 A. Average or composite spectra derived from large 
survey databases involving hundreds of AGN (LBQS: Francis et al. 1991; HST: Zheng et al. 1997; 
1998; FIRST: Brotherton et al. 2001; Sloan survey: Vanden Berk et al. 2001) make it possible to 
generate a "typical" quasar spectrum from blueward of Lya to Ha. 

What is unclear is whether such composite spectra are astrophysically useful beyond, per- 
haps, allowing us to identify lines that are too weak to be seen in individual source spectra. The 
fundamental question is whether the similarities or the differences in AGN line and continuum 
phenomenology tell us more about the underlying physics. The Eigenvector 1 concept has been 
advanced as a possible "H-R Diagram" for AGN in the sense that it appears to provide parameter 
space discrimination between all major classes of broad line sources as well as a correlation for, 
at least, radio quiet sources (Sulentic et al. 2000b). The correlation and distribution of radio-quiet 
sources (with FWHM(H/3)< 4000 km s -1 ) in El have been reasonably well fit with a model that 
sees accretion rate, convolved with source orientation, as the principal physical driver (Marziani 
et al. 2001). The generality of El, of course, requires much further testing but unless it is far 
off the mark, it suggests that average or composite spectra should be viewed in the same way that 
one would view an average stellar spectrum taken over the full effective temperature range that 
is observed in the H-R Diagram (spectral types O-M). Interpretation of composite spectra from 
heterogeneous samples will be complicated in two ways: (1) they will be subject to selection bi- 
ases dependent on the relative number of sources with each "spectra type" in a given sample and, 
more significantly, (2) they will average the spectra of sources with dramatically different phys- 
ical properties. We suggest that the most useful approach to averaging AGN spectra lies within 
the El context. We present average spectra for fixed domain quadrants in El followed by a brief 
discussion of important differences. The emphasis is on the FWHM(H/9) measures although large 
differences in the equivalent width parameter are also seen. Full discussion of the equivalent width 
differences must wait for significant numbers of very high S/N measures for Fell op t weak sources. 



2. Sample, Data Reduction and Analysis 

Our current sample includes optical spectra for 187 sources obtained between 1988 and 1996. 
They cover the wavelength range that includes H7, H/3, and most of the optical Fell opt emis- 
sion (4200-5700 A in the rest frame). Spectra were obtained with the following telescopes and 
spectrographs: ESO 1.5m (B&Ch), San Pedro Martir 2.2m (B&Ch), Calar Alto 2.2m (B&Ch), 
KPNO 2.2m (Gold), and Mount Ekar (Asiago) 1.82m (B&Ch). Data for 52 objects can be found 



-3- 



in Marziani et al. (1996) and the unpublished part of the dataset will appear in a forthcoming paper 
(Marziani et al. 2002). Spectra were taken with very similar instrumental setups including typically 
120 A/mm dispersion and 2 arcsec slit width yielding resolution in the range 4-7 A FWHM. The 
S/N is typically in the range ps 20-30. Spectra with S/N < 10 were excluded. The sample includes 
spectra for 187 type 1 AGN with redshifts z < 0.8, my < 16. Table 1 summarizes the properties 
of our data sample within the context of Eigenvector 1 (see §3.1 for bin definition). Composite 
spectra were obtained by taking the median of all spectra with lAbinning. The S/N of the com- 
posite spectra was > 100 except for bin A3 where only 5 objects were available, and S/N?«50. No 
composite spectra were obtained for the rare objects in the outlier domain (FWHM(H/3bc)> 4000 
km s _1 and RF e ii> 0.5) which spans several El bins. 

The following analysis procedures were applied to each spectrum (see Marziani et al. 1996 
for details). (1) Deredshifting using the narrow line component of Hf3, and/or [Om]A4959,5007 
as rest frame measures. (2) Normalization using the best estimate for a local continuum at A « 
5100 A. (3) Subtraction of Fell A4570 emission blends using a template based on the spectrum 
of IZwl which also yielded the Fell blend equivalent width W(Fell A4570). (4) Subtraction of 
the [Om]A4959,5007 and HellA4686 lines. (5) Subtraction of the narrow component from H/3 
(H/9 NC ) using a Gaussian model. In cases where no clear inflection between broad and narrow com- 
ponents was seen, two alternative approaches were employed: (a) we assumed that FWHM(H / 9 NC ) 
« FWHM([Om] A5007) or (b) we assumed that FWHM(H/3 BC ) ~ FWHM(Fell A4570). The latter 
approach was useful for Narrow Line Seyfert 1 sources. Following these procedures we extracted 
H/3 BC and derived equivalent width and FWHM measures which allowed us to populate El. Host 
galaxy contamination was not significant for the luminous AGN included in this sample. The 
typical uncertainty for FWHM(Hftc) is about 10% and for W(H/3 BC ) and W(Fell A4570) about 
10-15%. Errors on R Fen are estimated to be less than 0.2. 



3. Results 
3.1. Spectral Types 

The distribution of 187 sources in the optical El plane (FWHM(H/9bc) vs - RFeii) is shown 
in Figure 1 . We binned the optical parameter plane as shown in Figure 1 and as defined in Table 
1. The bins were arbitrarily set to constant A FWHM(H/3 B c) = 4000 km s _1 , and A R Fe ii= 0.5. 
A finer subdivision is not warranted at this time by the accuracy of the measures or the hetero- 
geneity of the data sample. The adopted A-B bin labeling reflects our earlier suggestion that two 
distinct radio-quiet AGN populations may exist: Pop. A ("pure" radio-quiet) and Pop. B (same 
as the radio-loud domain; see Table 1). The numeral accompanying the letter designation for A 
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bins reflects the increasing strength of PvFeii- The "+" designations for B bins reflect increasing 
FWHM(H/3 BC ). There are many measures that distinguish between radio-quiet Pop. A and B as 
defined with a boundary at FWHM(H/3 BC ) ~ 4000 km s _1 . There are also many measures that 
support a phenomenological commonality between population B and radio-loud sources (Sulentic 
et al. 2000a,b,c). Table 1 summarizes the properties of domain space bins that are occupied by a 
significant number of sources. The composite (median) spectra after processing step 2 (continuum 
normalized) and after step 3 (Fell opt subtracted) are presented in the left and right panels of Figure 
2 respectively. The median spectra are available in digital format on the World Wide Web. 6 

It is important to note that the El binned median spectra are unlikely to be strongly influenced 
by any luminosity dependence since: (a) there is, so far, no convincing evidence for a difference 
in line profile properties with source redshift/luminosity (Sulentic et al. 2000a) and (b) for any 
sub-sample, AM B < 1.5 < om b (see Table 1). It is also worth noting that the median M B of our 
sources is close to the limit separating Seyfert 1 galaxies and quasars (M B =-23.0). This should 
not be considered a particular problem because there is no discontinuity in the AGN luminosity 
distribution. Radio-loud sources in our sample show a higher mean luminosity perhaps due to the 
inclusion of many beamed sources but Pop. B radio-quiet sources, co-spatial with them in El, are 
not more luminous than Pop. A radio-quiet sources. 

Our W(Fell A4570) uncertainty estimate (§2) is valid if Fell A4570 is clearly detected above 
the noise. However, the minimum W(Fell A4570) for which it is possible to visually detect Fell 
A4570 depends both on S/N and on Fell A4570 FWHM. For S/N « 20. upper limits for Fell A4570 
detection is estimated to be m 14, 17, 20A for FWHM = 2500, 5000, 7500 km s" 1 respectively. 
The utility of the composite spectra is also evident because of the much higher S/N in bins Bl (and 
Bl + ), which contain many upper limits because Fell A4570 is intrinsically weak in Population B 
sources. In those bins it is possible to obtain more reliable W(Fell A4570) measurements from 
composite spectra. 



3.2. Median H(3 B c Profiles Across El 

Median spectra of the H/3 BC profile are shown in Fig. 3. They were extracted from the median 
spectra after Fe ll op t and continuum subtraction to take advantage of the much higher S/N of the 
compositer spectra. We derived a pure H/3 BC profile by first subtracting [Om]A4959,5007 and 
HellA4686 and then fitting H/3 BC with a high-order spline function to avoid the effects of finite 
S/N and to interpolate across the H/3 NC . This approach has proven to be successful for measuring 



6 At the address http://panoramix.pd. astro. it/^marziani/medians. zip 



-5- 



line parameters in noisier line profiles (Marziani et al. 1996). It is preferable to template fitting 
because it makes no assumptions about the shape of the profile. Table 1 presents FWHM(H / 9 B c) 
as well as line centroids at 1/2 and 1/4 maximum intensities derived from the median spectra. The 
median H/3 B c profiles in Pop. A are almost symmetric with a slight blueward asymmetry in bin 
A2. The blue asymmetry becomes visually apparent in bin A3 which contains sources that are 
extreme in many ways. Table 1 shows that the profiles of NLSy 1 do not differ from bin Al and A2 
median spectra. The smaller median FWHM H/3 B c reflects an arbitrary cutoff (<2000 km s _1 ) in 
the definition of NLSy Is. 

The profiles of Pop. A sources are very different from those of Pop. B. The median pro- 
files of population B sources are red asymmetric with the strongest asymmetry in the (few) bin 
Bl ++ sources. The different shape may point towards different BLR structure/kinematics in these 
sources. Profile shape was not selected a priori in the population definition that was based on a 
possible FWHM(H/?bc) break. An alternate way to describe profile asymmetry involves fitting a 
functional form to the median profiles. The H/5bc profile in NLSy 1 sources is well fit by a Lorentz 
function (see also Veron-Cetty, Veron, & Goncalves 2001). The same is generally true for Al and 
A2 median profiles while A3 H/^bc is based on only five sources. 

A Lorentzian function is definitely not a satisfactory description for bin B median profiles. 
The H/3bc profile is significantly different and this can be seen even in Fig. 2. The profile be- 
comes more Gaussian but we were unable to fit bin B profiles with any simple and physically 
meaningful function. The best fits to H/3 B c bin Bl and B1+ profiles required the sum of two 
Gaussians: (l)a FWHMm 4000-5000 km s~ 1 unshifted Gaussian core and (2) a broader FWHMm 
10000km s^redshifted (Av r ~ 1000 km s^ 1 ) Gaussian base. The main difference between Bl 
and Bl + may be related to the strength of the "very broad component". We will show in the next 
section: (1) that many radio-loud and some radio-quiet Pop. B sources show unambiguous very 
broad components and (2) that the interpretation of the profile in terms of two components may 
have a straightforward physical meaning. 



4. Discussion 

We find significant and systematic differences in median spectra across El. FWHM(H/9bc) 
and Rfcii both show systematic changes and the other El parameters reinforce this dichotomy 
(pop. A sources show a systematic blueshift of the CIVA1549 broad line profile and and a soft 
X-ray excess while pop. B source do not: Sulentic 2000ab). Pop. A sources show a more symmet- 
ric profile with the largest R Fen values (indicating large W(Fell A4570) and somewhat depressed 
W(H/?bc))- As one enters the population B domain the radio-loud source fraction increases. H/9bc 
becomes increasingly broad and red asymmetric for both radio-loud and radio-quiet sources. Over- 
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all NLSyl/Pop. A spectra show lower ionization spectra than Pop. B (Marziani et al. 2001). 
Comparison between high (CivA1549 B c : HIL) and low (H/? BC : LIL) ionization emission lines 
suggests that they are emitted in disjoint regions in pop. A. Contrarily, both LIL and HIL may 
arise from a single region in Pop. B sources (Marziani et al. 1996: Sulentic et al. 2000a,b). It is 
more accurate to say that we cannot rule out the possibility that both HIL and LIL arise from the 
same line emitting regions (broad and very broad) in population B (Sulentic 2000a). The H/^bc 
profile analysis indicates that the BLR - and specifically the region emitting the low ionization 
lines - is most probably structurally different in Population A and Population B sources. It is be- 
yond the scope of the present paper to report on detailed model calculation. We note however that 
the Lorentz profile is consistent with emission from an extended accretion disk. This reinforces 
the suggestion that the LIL spectra in Pop. A sources arise from a disk. The situation is less clear 
for Pop. B sources where the Eddington ratio may be much lower. There is good evidence that 
sometimes only one of the two emission components is present in Pop. B sources (pure BLR or 
pure VBLR; see Sulentic et al. 2000c). In a companion paper, we suggest that even the Narrow 
Line Region structure/kinematics may be changing along the El sequence (Zamanov et al. 2002). 



4.1. The Very Broad Line Region Issue 

Can the double Gaussian model that is needed to fit Pop. B (and radio-loud) profiles be phys- 
ically justified? Several lines of evidence point towards the existence of a Very Broad Line Region 
(VBLR) at the inner edge of the BLR (Corbin 1997a,b). Emission from this region may be thought 
of as a sort of inner, large covering factor, "boundary layer" where gas begins to become optically 
thick (Marziani & Sulentic 1993; Sulentic et al. 2000c). In objects like the luminous quasar PG 
1416-129, almost pure VBLR emission survived after a continuum intensity decrease effectively 
quenched the classical BLR. 3C232 (Marziani et al. 1996) may represent a good radio-loud analog 
of PG1416. Other examples of possible pure VBLR line sources can be found in Sulentic et al. 
2000c. HellA4686 has also been known for some time to show a profile systematically broader 
than H/9 in some sources (see Osterbrock & Shuder 1982) and this can be seen in Figure 2. Un- 
fortunately HellA4686 is usually weak and/or blended with stronger H/3 and Fell emission (it is 
on the red end of Fell A4570). Detection of the HellA4686 broad component is therefore a func- 
tion of Fellopt strength and S/N. The strongest very broad HellA4686 feature is seen in the bin Al 
median spectrum where we have a favorable combination of S/N, H/3 width, Fell opt weakness and 
HellA4686 strength. The latter bias reflects the inclusion of data from an observing run dedicated 
to studying strong HellA4686 emission. A very broad HellA4686 profile is sometimes seen in pop. 
A sources (e.g. Ton S 180) with much narrower Fell A4570However, in I Zw 1, the highest S/N 
spectra fail to detect any trace of HellA4686 emission. At this time the question of HellA4686 
emission must be considered on an object-by-object basis. We note that HellA4686 strength was 
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an orthogonal (eigenvector 2) parameter in the principal component analysis of Boroson & Green 
(1992). 

5. Conclusion 

We generated average QSO spectra in the 4200-5700 A spectral region for fixed optical pa- 
rameter bins in El . In fixed parameter range bins the mean spectral type appears to change system- 
atically across El. This systematic behaviour is reflected in El X-ray and UV measures as well. 
Composite spectra in an El context should provide useful input for theoretical modeling. Ignoring 
the diversity of the AGN H/3 B c profile may make any theoretical modeling of the BLR kinematics 
and dynamics unreliable. 

The authors acknowledge support from the Italian Ministry of University and Scientific and 
Technological Research (MURST) through grant and Cofin 00-02-004. 
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Table 1. Sample Properties 



Type 


RfoII 


FWHM(H/3 B c) 
[kins- 1 ] 










FWHM(H/3 BC ) b 
[kms- 1 ] 


R-FcII 


c(l/2) b - c 
[kms- 1 ] 


c(l/4) b - c 
[kms" 1 ] 


Best Fit 


Total 


all 


all 


187 


64 


-23.7 


2.0 




0.37 








NLSyl 


all 


<2000 


24 


2 


-23.0 


1.7 


1500 


0.61 


40 





Lorentz 


Pop A d 


all 


<4000 


80 


11 


-23.3 


1.8 


2300 


0.49 


20 


30 


Lorentz 


A3 


1.0-1.5 


- 4000 


5 





-24.5 


2.2 


2350 


1.23 


-280 


-290 


Lorentz 


A2 


0.5-1.0 


- 4000 


26 


3 


23.3 


1.5 


1950 


0.70 


-20 


-30 


Lorentz 


Al 


0.0 - 0.5 


- 4000 


49 


6 


-23.2 


1.8 


2250 


0.31 


30 


250 


Lorentz 


PopB c 


all 


>4000 


97 


50 


-24.2 


2.0 


5700 


0.22 


200 


700 


Double Gaussian 


Bl 


0.0 - 0.5 


4000 - 8000 


76 


35 


-24.2 


2.0 


5600 


0.23 


250 


700 


Double Gaussian 


B1 + 


0.0 - 0.5 


8000- 12000 


17 


12 


-24.1 


2.1 


10000 


0.17 


270 


1700 


Double Gaussian 


B1++ 


0.0 - 0.5 


> 12000 


4 


3 
















outliers 


>0.5 


>4000 


10 


3 


-23.0 


2.3 




0.80 









a Absolute B magnitude from the Veron-Cetty and Veron Catalogue, for Ho = 50 km s -1 Mpc -1 , go = 0. 
b Values measured on a high order spline function representative of the H/3 B c of each median spectrum. 

c Centroid at fractional intensity | are defined as c(i/4) = (Ag + A/j — 2Ao)/Ao, where Ao = 4861.33 A, and As and A^j are the wavelength of the blue 
and red line side at the given fractional intensity. 

d Pop. A = Al U A2 |J A3. 

e Pop. B =B1|JB1+ U B1++ » with outliers (O) excluded 
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Fig. 1. — The optical parameter plane (Rf c ii-FWHM(H / 9bc)) of Eigenvector 1. This is the largest 
sample that yet displayed in an El context. Dotted lines correspond to the borders between dif- 
ferent spectral types also listed in Table 1. The abcissa is R FcII =W(Fell A4570)/W(H/3 BC ). The 
thin solid line indicates the boundary of "outlier" sources and is, so far, a zone of avoidance in 
quasar parameter space. Horizontal arrows denote Fell A4570 upper limits. Note that upper limits 
dominate single Rp e n measurements in the Bl and Bl + domains, while the higher S/N ratio of the 
median spectra allow for a more accurate determination. 
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Fig. 2. — Left panel: Composite (median) spectra of the El parameter bins defined in Table 1 and 
Figure 1 . Right panel: the same composite spectra with Fell op t emission subtracted. 
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Fig. 3. — Fellopt, and continuum subtracted H/3 composite line profiles for the different El pa- 
rameter bins defined in Figure 1. Solid lines show the profile H/9bc after subtraction of narrow 
line emission components. A Lorenztian fit (dashed line) is superimposed on the NLS1, Al, A2 
profiles. The individual components of a double Gaussian (thin lines) and resultant fit (dashed line) 
are shown for Bl and Bl + . 



